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We consider the dispersion energy of a pair of dipoles embedded in a metallic waveguide with 
transverse dimension a smaller than the characteristic dipolar wavelength. We find that a sets the 
scale that separates retarded, Casimir-Polder-like, from quasistatic, van der Waals-like, interactions. 
Whereas in the retarded regime, the energy decays exponentially with inter-dipolar distance, typical 
of evanescent waves, in the van der Walls regime, the known free-space result is obtained. This short- 
range scaling implies that the additivity of the dispersion interactions inside a waveguide extends 
to denser media, along with modifications to related Casimir effects in such structures. 
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I. INTRODUCTION 

Electric dipoles interact through their electromagnetic 
fields. Hence, the propagation and scattering proper- 
ties of these fields in the medium where the dipoles are 
embedded, strongly influence the spatial dependence of 
dipole-dipole interactions. In the framework of quantum 
electrodynamics (QED), these interactions are described 
as being mediated by virtual photons in different field 
modes. The space dependence is then determined by the 
spatial wavefunction of the field modes and hence on the 
geometry. 

When one of the dipoles is initially excited, the result- 
ing interaction induced by virtual photons is the so-called 
resonant dipole-dipole interaction (RDDI) , where the ex- 
citation is dynamically exchanged between the two iden- 
tical dipoles [TJ[2]. In free space, where the photon modes 
are plane waves, RDDI scales as 1/r 3 for inter-dipolar dis- 
tances r, much shorter then the typical dipolar-transition 
wavelength A, and cos(27rr/A)/r for r ^> A [T]. Stud- 
ies of other geometries however, where RDDI is medi- 
ated by confined photon modes, e.g. one dimensional 
surface plasmon polaritons [3J, photonic band gaps 0], 
cavity modes [S] and hollow metallic waveguides [7] to 
name a few, yielded rather different space-dependencies. 
For example, in a rectangular hollow metallic waveguide 
(MWG), depicted in Fig. 1, long range RDDI at dis- 
tances of 100A was found possible, by exploiting the ex- 
istence of cutoffs in the photon spectrum [6]. For modes 
with cutoff frequencies higher than 2ttc/X, the interac- 
tion is mediated by evanescent modes, i.e. decaying ex- 
ponents, that sum up to give the free-space result 1/r 3 
at distances much shorter than the size of the waveguide 
confinement (a and b in Fig. 1) [7J. 

Consider now the van der Waals (vdW) dispersion in- 
teraction, namely when both dipoles are in their ground 
states, and their mutual interaction is mediated by vir- 
tual photons from the vacuum. In free space and for short 
distances w.r.t the dipole wavelength, r <C A, the vdW in- 
teraction scales as 1/r 6 , a direct consequence of the qua- 
sistatic limit of the RDDI, 1/r 3 . However, at longer dis- 
tances, r 3> A, in the so-called retarded regime, Casimir 
and Polder obtained a 1/r 7 dependence 2.8,9 . Different 



spatial scalings were found however, when two interacting 
bulk objects rather than point dipoles were considered. 
Casimir studied the case of two parallel metal plates and 
found an attraction force that scales as l/d 4 , with d the 
distance between the plates [TO] . Since then, related ef- 
fects are constantly being studied, mainly considering the 
interaction between dielectric or metallic objects of var- 
ious geometries [3 [TT]. Such experimental [T2"lu"5] and 
theoretical [T6H20] studies demonstrate the role of retar- 
dation and geometry dependence of these vdW-related 
phenomena |21j . 

Here however, we would like to take a somewhat differ- 
ent point of view towards the spatial-dependence of the 
dispersion interaction, reminiscent of that we described 
for RDDI. Namely, instead of considering objects with 
different geometries and calculate their interaction en- 
ergy, we refer back to the original vdW configuration of 
a pair of point dipoles while changing the geometry of 
their surrounding environment, and hence the structure 
of the mediating virtual photon modes. Specifically, here 
we study dispersion interactions between two dipoles in- 
side a metallic waveguide, in the case where the typical 
dipole wavelength A is much larger than the typical trans- 
verse confinement of the waveguide (a and b). 

Our results, obtained analytically by QED perturba- 
tion theory, are analogous to those of RDDI in a MWG 
O [7J . They are rather general although full expressions 
are then derived for the case of a rectangular MWG. Since 
the smallest cutoff frequency for the mediating transverse 
modes of the waveguide is of order ctt /a, they become 
evanescent for the dipole frequency 2nc/X < en /a. In- 
deed, we obtain the dispersion energy as a sum of de- 
caying exponentials, J2 n bne~ KnZ , where z is the inter- 
dipolar distance on the propagation axis of the MWG, n 
is an index of a transverse mode, and b n , K n are constants. 
Two distinct regimes are than recognized: the retarded, 
wave-like, limit, at long distances z 3> a, where the in- 
teraction decays like e~ K ° z , with ko the smallest of all 
the k„'s, and the vdW, short-distance limit at z « o, 
where the free-space result 1/r 6 is restored. Hence, as 
opposed to the free-space case, where A sets the scale 
for retarded (Casimir-Polder) and non-retarded (vdW) 
regimes, here we find that the relevant scale is a. Two 
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FIG. 1: Dispersion interaction inside a rectangular hollow 
metallic waveguide. The transverse x, y dimensions of the 
waveguide are a and 6, respectively. A pair of interacting 
dipoles, represented by black dots, are located inside the 
waveguide and along its propagation axis z. 



consequences then emerge: (1) unlike free-space, retar- 
dation effects are dominant here at distances set by a, 
much shorter than the wavelength A, and (2) the inter- 
action has a range much shorter than in free-space due to 
the exponential decay in the retarded regime. The latter 
also suggests that the approximate additivity of dipolc- 
dipole and dispersion interactions in a dilute system of 
many interacting dipoles [9], is expected to hold here at 
even higher densities. 

Another interesting aspect of our results is related 
to the character of the dominant transverse modes in 
the Casimir-Polder and vdW regimes. Whereas, the 
smallest cutoff frequency is obtained for the lowest or- 
der transverse-electric (TE) mode and hence it is domi- 
nant in the retarded limit, the transverse magnetic (TM) 
modes are dominant in the short-distance limit and their 
contributions sum up to give the familiar free-space re- 
sult. 

The paper is organized as follows. In Section II we 
provide a general discussion and formalism for the cal- 
culation of the dispersion energy mediated by photon 
modes that posses a cutoff, and introduce the modes of 
the MWG. Next, in Section III, we present the calculation 
and results for the dispersion energy in a MWG, whereas 
the comparison with free-space in the corresponding vdW 
(short-range) and Casimir (retarded) regimes is analyzed 
in Section IV. Our conclusions are finally given in Section 
V. 



II. DISPERSION INTERACTION VIA MODES 
WITH CUTOFFS 

A. The system 

We consider two identical dipoles, e.g. made of atoms 
or molecules, with ground state \g) and excited levels 
{|e)} with corresponding energies E e w.r.t that of the 
ground state. The excited levels possess dipole-transition 



matrix elements d e to the ground state. The dipoles are 
placed inside a hollow waveguide that supports photon 
modes with normalized mode functions, 

<nk^) = ^ ik ^ mnk (x,y). (1) 

Here k is the wavenumber at the propagation axis z with 
quantization length L, fi is the index of the polarization 
and m and n are indices of the transverse mode whose 
dependence on the transverse coordinates is described 
by fe (a;, y). Each of the transverse modes ix mn has a 
dispersion relation 

= c V k ?nn + k2 , ( 2 ) 

where c is the speed of light in the material that fills the 
hollow waveguide and ck mn is the cutoff frequency of this 
mode. The Hamiltonian of the two dipoles 1,2 and the 
field modes is then described by 

Ha = ^2 E M(e v \, 

f=l,2 e„ 

U — \ ^ hi if 

12 F — " ,UJ mnk u mnk u mnk> 
^mnk 

2 

H AF = -HJ2J2 E (l^)(5|+h.c.)(i ff ™ fe a^ nfe +h.c), 

(3) 

where = d e , y • < nnfc (r l/ ) is the dipole cou- 

pling between the |e) -s-J- \g) transition of atom v, located 
at r„, and the /x mn fc mode with lowering operator <z^ nfc , 
and e being the permittivity inside the MWG. 

B. QED perturbation theory 

In order to calculate the dispersion interaction energy 
we follow the QED perturbative approach of Ref. [2]. 
The interaction energy U is obtained as the fourth order 
correction to the energy of the ground state of the system, 
\G) = 1 51,52,0) with energy Eq = 0, where both atoms 
are in their ground state and the photon modes are in 
the vacuum |0), 

(E h - E G )(E l2 - E G ){E l3 - E G ) 

(4) 

Here \Ij) are intermediate (virtual) states, and E q is the 
free Hamiltonian (Ha + Hp) energy of the state \q). The 
above sum contains 12 possible terms, each of which de- 
scribes a different virtual process and can be represented 
by a diagram. In Fig. 2 we present the diagrams of two of 
these processes: the one in (a) includes the intermediate 
states |ix) = |ei,3 2 , l„„ fc ), \h) = |ei,e 2 ,0) and |7 3 ) = 

lfi' e 2> Im'n'k')' with \ l mnk) = a m„fcl°> > whereas the one 
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FIG. 2: Dispersion energy calculation by QED perturba- 
tion theory. Two of 12 possible processes that contribute 
to the energy correction of the state \G) = Iffi, <?2, 0), Eq. 
Q. (a) Dia gram of the process that includes the inter- 
mediate states = |ei, l^ nfe ), |Ja) = |ei,e 2 ,0) and 

|/a) = Isi.ea.l^), with = a£ nfc |0>. Its corre- 

sponding term in the sum Eq. p} has the denominator 
D a from Eq. |5|. (b) Here the intermediate states are 

\h) = |ei,ff 2 , C k ), I/2) = |fli,fl2,l{Ufcl 
Iffli e2, lm',j'fc'}> with the denominator D;, from (|5 



mnfc^m'n'fe'L aIld l J 3> = 



in (b) describes the process — \e\,g<z,^ 
l5i,32,C nfc 1 m'„'fe'> and l J 3> = Iffi.ea,!" 



m 11k I 



L mnk J ~m' n' k' I a,1AVX I 3 

to show that all 12 terms in the sum 



n'k')- It is easy 
have the same 



numerator, but are distinct in their denominators. For 
example, the terms that correspond to the diagrams in 
Figs. 2(a) and 2(b) possess the denominators, 



Fig. 2. In fact, there are three more diagrams with de- 
nominator D a and in a similar fashion one can show that 
D a is smaller than the rest of the possible denominators 
of the remaining 7 diagrams, out of the total 12. Hence, 
the dispersion energy in the tight-confinement case can 
be approximated by considering only the contribution of 
the diagrams with D a . Namely, we obtain the energy U 
by 4 times the contribution of diagram (a), 



U 



^e.Q_i de i j d e 2 1 d e 1 q 



' (2^e) 2 EE' £ , + E ' 

ei,e 2 i]lq 1 i 



F e2 ijF ei lq, 



F,, 



mn,eij 



mn,eij ' 



fimn 
f oo 



dk 



^mnk 



X E mnk,i (^2 , V2)E^ nk] {x llVl )e 



ikz 



(7) 



Here z is the inter-dipolar distance in the waveguide 
propagation axis z, whereas (x u ,y v ) is the transverse po- 
sition of dipole v. The indices i,j,l,q run over the pro- 
jections of the vectors d e and E^ nfe (a;, y) onto the x, y, z 
directions, e.g. d ex — d e ■ e x . Further simplification can 
be made by applying the inequality |6]) in the above ex- 
pression for F^ n eij , obtaining 



mn,ezj 



dkE L l kA x 2,y2)E^ nk Ax 1 ,y 1 y 



ikz 



(8) 



D a = {hu» 

(tkj, 



Dh 



mnk 



respectively. 



E ei )(E ei - 



fi^m'n'k 



m' 'n'k 



Ee 2 ), 



Eq. ([7]) show that the spatial dependence of the disper- 
ses )i sion energy U is encoded in F^ nei -, while Eq. |8j) fur- 
(5) ther simplifies its calculation by the Fourier transform 
of a multiplication of the waveguide modes' transverse 
profiles. 



1. The case of tight confinement: a,b <^ 2irhc/E e 

Let us now consider the case where the typical wave- 
lengths associated with all the dipolar transitions, A e = 
2irhc/E e is much larger than the transverse confinement 
of the waveguide, a and b. Since the cutoff wavenumbers 
of the transverse modes, k mn , result from the confine- 
ment, they must satisfy k mn > n/a or ir/b V/i m „, and 
we obtain in our case k mn 3> E e /(hc) for all states |e) 
and modes \x mTl . Then, considering also the dispersion 
relation, Eq. §2§, wc find 

a; m«fe >£ 'e/^ Ve,n mn ,k. (6) 

This means that D a <C Di, such that the contribution of 
diagram (b) is negligible w.r.t that of diagram (a) from 



C. The modes of a MWG 



We consider the rectangular hollow MWG shown in 
Fig. 1, whereas any of our calculations can be easily gen- 
eralized to e.g. a circular hollow MWG. We also take the 
limit of a perfect metal, which is reasonable when the rel- 
evant frequencies E e /h, cir/a and cir/b are small enough 
w.r.t the plasma frequency of the metal. There are two 
possible polarizations for the transverse modes, TE and 
TM [25], and by normalizing their transverse profiles, 

i.e. demanding f£ dx J Q & dy^ nk ■ (x,y)^ nk {x,y) = 1, 
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(9) 



where A = ab is the transverse area of the waveguide, and 
we note that the transverse profile of the TE modes, E^, 
does not depend on the wavenumber k (or on frequency), 
hence the index k is omitted. The mode frequencies a/* ■ 
with /j = TM, TE are those from Eq. (|2) with the cutoff 
wavenumbers k mn given by 



k mn = y/ {mil /a) 2 + (nix/b) 2 



(10) 



The indices to, n are positive integers, whereas the lowest 
modes are TM U and TE Q1 or TE W for TE and TM, 
respectively. 



III. ENERGY CALCULATION AND GENERAL 
RESULTS 



In this section we present the calculation and general 
results for F^ ln ei - form Eqs. (|7j) and |8j), for all dipole 
orientations i, j and mode polarizations (i = TM, TE. 



which cannot exceed k^l — \J a 2 + b 2 jix. The exponen- 
tial decay is a consequence of interaction via evanescent 
modes, since all of the frequencies of the dipole transi- 
tions E e /h, were assumed here to be lower than all the 
cutoff frequencies of the MWG, ck mn (see Sec. IIB). 

Turning to the i, j = x, x case, we insert Ej^ k x (x, y) 
from Eq. ([9| into Eq. and obtain, 



TM 4 / TO7T 

^mn,exx = 7™t»l I T I cos 



x / du- 



^UK mn Z 



(mix \ 2 (nix \ 

{-TV sm {TV 

(13) 



u 2 + 1 



We perform the integral by contour integration with reg- 
ularization and get, 



j?™ - iZU f mn 

J± \ /».,,, ,, fX 



cos ^ a; j sm ) e 



(14) 

obtaining again exponential decay with the exponent 
kmn- We proceed in a similar fashion for the rest of the 
options and find 
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2mix \ . / 2mr \ ., 
— x sin — y )e - 
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and F„ 



TM _ rpTM 



F mn,en for (^l.J/l) = (^2,2/2) = 



(15) 



TM modes 



We begin with the i, j — z,z case, namely insert the z 
component of the TM mode E™ hz (x, y) = E™ k (x, y) ■ 
e z from Eq. ([9| into the integral in Eq. ([8|. We get, 



F 



TM 



mn.ezz 



A 



2 (mix \ 2 (nix \ 



1 



(11) 



TE modes 



By virtue of Eq. (|9j) we note that E^^(x, y) (i — x, y) 
is independent of k , such that it can be taken out of the 
integral in Eq. (7| for F^ mj . Then, by using Eq. 
we are left to solve the integral, 



TTTi 



1=1 du- _ 

-oo V" 2 + l + U e 



(16) 



where u is a dimensionless integration variable. Here 
we assumed for simplicity (xi,y±) = (x2,V2) = ix,y). 
The above integral is simply the Fourier transform of 
a Lorentzian, giving rise to an exponential decay with 
distance z, 



F, 



TM 



4tx 

T 



. 2 / mix \ . 2 (nix \ 
fc mn sm y— — xjsin ^-^-J/J 



,-fe« 



_ (12) 

The typical range of the TA/ m „-mediated interaction be- 
tween the z component of the dipoles is than set by 



with C = k mn z and u e = E e /(hck mn ). By contour inte- 
gration methods, this integral is found to be equal to 



I = ~2u e 



du 



VvTT 



l + ul 



/i TT^i 



-(u _ 



-2u e K (()(17) 



where K Q (u) is the zeroth order modified Bessel function 
and where the lowest order approximation due to the 
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tight confinement, u e <C 1, was taken. We thus obtain 
for any i = x,y and j = x, y, 



F mn,ei 3 = E 1 m ^ l {x 2 ,y 2 )E^ J {x 1 ,y 1 ) 1 ^K (k mn z). 



-jTE 



E„ 



he 



(18) 



IV. 



COMPARISON TO FREE-SPACE: VDW 
AND CASIMIR-POLDER REGIMES 



After we obtained the general expression for the dis- 
persion energy for every polarization of the dipoles and 
mediated by any of the transverse modes fx mn , we now 
wish to consider two relevant limits; i.e., the retarded, 
Casimir-Polder limit, and the quasistatic, vdW limit, and 
compare the results with their free-space counterparts. 



A. Casimir-Polder: z 3> a,b 

In the tight confinement case we consider, the scale for 
the on-axis distance z at which wave effects emerge is set 
by the confinement a, b. However, since a, b <C 2irhc/ E e , 
wave effects in the confined case mean decay rather than 
oscillation with z due to the evanescent nature of the 
modes. 



TM contribution 



Eqs. (12 1,(14) and (15) show that for any dipole ori- 



entation, the TM-mediated energy decays exponentially 
with the exponent k mn like, 



exp(— 7r- y/r 



(19) 



where b ~ a is assumed here for simplicity. Clearly, when 
z>a, the contribution of the lowest order TM mode, 
namely m = 1 and n = 1, overwhelms that of the rest 
of the TM modes, and the \x = TM contribution to the 
sum for F ei j in Eq. Q, can be taken solely as 
where 



TM 

eiji 



TP 1 
^11 



TM 



oc e 



a 



(20) 



2. TE contribution 



The argument of the modified Bessel function in Eq. 



(18), is exactly that of the exponent in (19), and in the 



retarded regime z 3> a, it is clearly much larger than 1 
for any m,n / 0,0 (as in TE). We then use the approx- 
imation Ko(k mn z) Ri •\/ 7r / (2&r) 



iz)e 



for k mn z 1. 



The contribution of the TE modes can then be approxi- 
mated by that of the two dominant lowest order modes; 
i.e. TEqi and TEiq, where 



7 TE 
10, eij 



e -Trz/a 

\J~nzja 



TE 



TP' — rsr 

) r 01,eij tx 



e -irz/b 
y/%z/b 



(21) 



3. Total dispersion energy 



From Eqs. (20) and (21 1 we conclude that the TE 



contribution is more dominant than the TM one due to 
its slower decay with z/a. The total energy U from |7|, 
in the retarded limit, can thus be approximated by that 
of the TEqx and TE W modes, 



U ~ (27re) 2 51 

ei,e2 ijlq 
i JpTE J?TE 



E ei + E e2 



l r 10, e 2 ij r 10, etlq 



;E 



TE 



P TE 

r 01,e 2 ij A 10,ei/ij 



P TE 
b 01 



<■/).) 



p TE 1 



(22) 



Using Eq. ([9]), we find 

2 7T 

E oi E i = ~^ sin ^ y ) 5i; 



E w E i = 7 =sin(-x)S ly , 

V A ci 

(23) 

where Sij is the Kronecker delta. In order to illustrate 
our results, let us take x\ = x 2 — yi = y% = x and a — b, 
and assume the dipoles are randomly oriented in space 
in a uniform distribution, i.e. (dg X ) = (d 2 ey ) — (<ig Z ) = 
(l/3)|d e | 2 . We then obtain 



U 



8tt^ 
9 



sin 4 fa) | 



v |d ei | 2 |d e2 | 2 



1 



A ei A e2 a 3 z 



-2irz/ a 



(24) 



where A e = hc/E e is the wavelength associated with the 
\g) f-> |e) transition. We note that the effect of the trans- 
verse position of both dipoles, x, is contained in the nu- 
merical prefactor sin 4 (^x) and becomes equal to 1 when 
the dipoles are in the center, x — a/2. 



4- Comparison to free-space result 

It is now instructive to compare the above result to its 
free-space counterpart, Uf s , in two different regimes. For 
a « z « A e Ve, the corresponding result in free-space is 
that of the non-retarded, vdW, limit [2] 
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24tt 2 e 2 E P 
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(25) 



where r, the inter-dipolar distance, is equal to z in our 
case. For a single excited level e\ — e 2 — e and x — 
a/2, the ratio between the MWG and free-space results 
becomes, 
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Ufs 



64tt 4 
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3 A? a 3 



(26) 



While the factor z 5 /(A 2 a 3 ) may become larger than 1, 
since z is between a and A e , the exponential decay still 
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Dispersion energy in the metal waveguide (MWG) U 
compared to its free-space counterpart Uf s , in the retarded 
z>o regime. Here a single dipole transition with wavelength 
X e is assumed, (a) a <C z <C A e : in this limit the free-space 
energy is in the vdW regime, Eq. 
from Eq. 



(251, and the ratio U/Uf s 
(b) 



100a. 



> A e 



( 26 I is plotted for \ c 
this li mit the free-space energy takes the C asim ir-Polder form 
[Eq. p7| ] and the ratio U/U fs from Eq. {26} is plotted for 
Ae = 10a. Both plots are linear in a semi-log scale, hence they 
show that the free-space energy is exponentially and many 
orders of magnitude larger than the short-range interaction 
med iated by the MWG modes, as suggested by Eqs. ( 26 1 and 
(28l. 



makes this ratio much smaller than 1, such that the free- 
space energy is much stronger than the one mediated by 
the MWG in this z limit. This is plotted in Fig. 3(a) for 
X e /a = 100. 

Next, we consider the limit z 3> A e Ve, where the free- 
space result takes the retarded, Casimir-Polder form [5] 
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(27) 



Taking again r = z, e% — e2 — e and x = a/2, the MWG 
to free-space ratio becomes, 
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1287T 
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(28) 



Although z 6 /(Aga 3 ) 3> 1, it is again the exponential de- 
cay that makes this ratio go practically to zero, as can 
be seen on Fig 3(b), for X e /a = 10. 

Eqs. (26 1 and (28) imply that the dispersion inter- 



action in a MWG is much shorter ranged than that in 
free-space. This is further discussed in the conclusions in 
section V. 



B. vdW: z < a,b 



not "sense" its metallic plates, free-space behavior is ex- 
pected to be observed. 



TM contribution 



The TM mn terms [Eqs . (|T2| ,([T4| and ([15])] possess 
an exponential decay e -^ m2 + n2 ^ z / a ( which in the z <C 
a, b is slowly decaying with m and n. Hence unlike the 
retarded regime, where only the TMu was considered, 
here we must consider all terms if we wish to account for 
arbitrary small z. In Sec. IVB3 below, we show how the 
summation over all the terms may yield the free-space 
result. 



2. TE contribution 



The z-dependence of the TE modes [Eq. (18)], goes 
like Ko(k mn z). It is obvious that for small enough to 
and n and in the limit z -C a,b discussed here, we have 
kmnZ € 1. In order to make our discussion clearer, let 



b, then k„ 



us assume a — u , mm n, mn 
to, n <C m* = a/(zn), we indeed have k 
take the approximation 



n-\/m 2 + n 2 , and for 

a v 1 

,z<l and can 



K (k mn z) w - ln(z/a) - ln(fc„ m a) + ln(2) - 7, (29) 

where 7 ss 0.577 is Euler's constant. 

Let us now try to estimate the total contribution of the 
TE modes: as z/a —> 0, the contribution of each to, n <C 
to* diverges like ln(z/a). Summing them together, we 
multiply ln(z/a) by roughly to* oc a/z terms, ending 
up with (1/z) In 2 divergence. If we wish to account for 
any, arbitrary small z/a, then to* — > 00 such that our 
estimation in fact includes the contribution of all terms. 
To conclude, we estimate 



/ j mn.eij 



lnz 

Z 



(30) 



Using this result in Eq. ([7]), we obtain that the vdW 
energy mediated by solely TE modes scales as U ~ 
(ln 2 z)/z 2 . Such divergence is much weaker then that 
of free-space, 1/z 6 , which we expect to obtain when the 
dipoles are very close. This implies, that the contribution 
of the TE modes to the dispersion energy at short dis- 
tances is not the dominant ingredient of the total energy 
in this regime. 



We define the vdW regime as the range of inter-dipolar 
distance z where wave effects are not significant. This 
is the case for z -C a, b, where the dipoles are close 
enough such that the exponential decay of the TE and 
TM evanescent waves, set by a and b, is almost not felt. 
Moreover, when the dipoles are placed in the center of 
the waveguide and are close enough such that they do 



3. Comparison to free-space result 

The above conclusion drawn for the TE modes contri- 
bution for the dispersion energy suggests that the TM 
modes are the dominant interaction mediators in the 
vdW, short-range, regime. We now reside to the cal- 
culation of the vdW energy due to TM modes. In order 
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FIG. 4: Evaluation of the sum in Eq. | [3l| : log-log plot of 
the result of the integral approximation, Eq. (321, compared 



to direct numerical summation. Excellent agreement is ob- 
served. 



to also relate it to its free-space counterpart, we assume 
that the dipoles are in the center, namely x = a/2 and 
y = 6/2, and take a = b for simplicity. Then, for the 
i,j = z,z component [Eq. (12)], we obtain, 



F,. 



V F™ 



An 2 



inn 

oo 



-F zz 



^ = E ^ 

m,n— 1 

w — y/ rn 2 -\-n 2 7T z I a 

x fe- 



rn 2 + n 2 sin 2 ( m- 



(31) 



We note that sin 2 (m7r/2) equals 1 for odd m and 
for even to, hence the sum F zz becomes F zz = 

For small z we 



^-jm^n—odd 

can approximate this sum by an integral, replacing 
J2 n =odd ~~ * (1/2) dn for n and to, and obtain 



4tt 2 ' 



1 

4n 2 (z/a) 3 ' 



(32) 



The integral approximation is tested by comparing it to 
a direct summation yielding excellent agreement, as can 
be seen in Fig. 4. 

In a similar fashion we calculate all other cases from 



Eqs. (12 1,(14 1 and (15 1 and get 



F 



F — F 

± EXT. ± < 



eyy 



1 

'2^3' 



F — F — F — 

J - exy — L exz — ± eyz u - 



(33) 



Inserting these results into Eq. ffi, we obtain the dis- 



persion energy in the vdW regime, 



U 



47r 2 e 2 z 



,6 1] Z! 



ei ,e2 ijlq 



E,. 



E e 



J i3,yy) 



3lq,zz 



2 ( S h,x 



(34) 



This has to be compared to the vdW result in free-space 
before averaging over dipole orientations, namely [2], 



U 



47r 2 e 2 



r 6 E E 



e l)^2 ijlq 



de2ideijde 2 ldeiq 

E ei + E e2 



°lq 



SRlRq 



i.e. R z = 1 and R x = R v 



(35) 



0. It 



where Ri = e,; ■ 
is easy to verify that Eqs. (|34| and ( |35[ ) yield identical 
results: for i,j,l,q — z the terms in the second row of 
both equations give 1, for i,j,l,q = x they give 1/4, for 
i,j = x and l,q — z they give -1/2 (repulsion) and for 
i 7^ j or I 7^ q they become 0. 

We conclude that in the vdW regime, namely z a, 
where no wave effects arc apparent, and when the dipoles 
are in the center of the waveguide, the TM modes are 
dominant and the dispersion energy is that of free-space, 
as can be expected. 



V. CONCLUSIONS 

In this paper, we applied ideas drawn from the modi- 
fication of resonant dipole-dipole interaction in confined 
geometries, to the study of dispersion interactions be- 
tween a pair of point-like dipoles. We considered a 
tightly confined waveguide structure whose transverse- 
modes cutoff-frequencies are all above the typical fre- 
quency of the interacting dipoles. We obtained analytical 
expressions for the dispersion energy at all inter-dipolar 
distances z. We found that the difference between re- 
tarded, wave-like (Casimir- Polder) , behavior and qua- 
sistatic (vdW) one, is set by the confinement length-scale 
a. For z <C a the interaction is mediated by TM modes 
and can become identical to its free-space counterpart, 
whereas in the retarded regime, z»a, the interaction is 
exponentially decaying with distance and is carried by 
TEqi and TEio evanescent fields. Hence, the result- 
ing interaction is much shorter range than that in free 
space. This has an implication on the non-additivity of 
the vdW and dipole-dipole interactions: in free-space the 
interaction energy between multiple dipoles is additive, 
namely, can be obtained by pairwise summation, as long 
as a/r 3 is small, where a and r are typical polarizability 
and inter-dipolar distance respectively [S]. Here however, 
this 1/r 3 scaling is expected to change such that even for 
higher densities, r < a 1 / 3 , the interaction may still be 
additive. This entails modifications of, e.g. the effec- 
tive dielectric constant of a gas of such dipoles, and all 



Casimir-associated phenomena in a MWG environment, 
which are based on the retardation and non-additivity 
[IT], that are both being altered here. 

Finally, let us address the generality of our results. 
Although they were explicitly derived for a rectangular 
hollow MWG, the same dependence on distance z is ex- 
pected for dispersion interactions mediated by any trans- 
verse modes with cutoff frequencies higher than that of 
the dipole transition. These may include, e.g. other 
cylindrical hollow MWG in various transverse shapes, 



tightly confined hollow-core fibers and high-order modes 
of electric transmission lines. 



Acknowledgments 

We appreciate useful discussions with Grzegorz Lach 
and Steven Johnson. The support of ISF and DIP is 
acknowledged. 



[1] R. H. Lehmberg, Phys. Rev. A 2, 883 (1970). 

[2] D. P. Craig and T. Thirunamachandran, Molecular 

Quantum Electrodynamics (Academic, London, 1984). 
[3] A. Gonzalez- Tudela, D. Martin-Cano, E. Moreno, L. 

Martin-Moreno, C. Tejedor and F. J. Garcia- Vidal, Phys. 

Rev. Lett 106, 020501 (2011). 
[4] G. Kurizki, Phys. Rev. A 42, 2915 (1990). 
[5] T. Kobayashi, Q. Zheng and T. Sekiguchi, Phys. Rev. A 

52, 2835 (1995). 
[6] E. Shahmoon and G. Kurizki, [arXiv:1205. 3064] (2012). 
[7] G. I. Kweon and N. M. Lawandy, J. Mod. Opt. 41, 311 

(1994). 

[8] H. B. G. Casimir. and D. Polder, Phys. Rev. 73, 360 
(1948). 

[9] P. W. Milonni, The Quantum Vacuum: An Introduction 
to Quantum Electrodynamics (Academic, 1993). 
[10] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 
(1948). 

[11] A. W. Rodriguez, F. Capasso and S. G. Johnson, Nature 

Photon. 5, 211 (2011). 
[12] A. O. Sushkov, W. J Kim, D. A. R. Dalvit and S. K. 

Lamoreaux, Nature Phys. 7, 230 (2011). 
[13] D. E. Krause, R. S. Decca, D. Lopez and E. Fischbach, 

Phys. Rev. Lett. 98, 050403 (2007). 
[14] C. C. Chang, A. A. Banishev, G. L. Klimchitskaya, V. M. 



Mostepanenko and U. Mohideen, Phys. Rev. Lett. 107, 
090403 (2011). 

[15] J. Munday, F. Capasso, V. A. Parsegian, Nature 457, 
170 (2009). 

[16] D. A. R. Dalvit, F. C. Lombardo, F. D. Mazzitelli and 

R. Onofrio, Phys. Rev. A 74, 020101(R) (2006). 
[17] T. Emig, R. L. Jaffe, M. Kardar, and A. Scardicchio, 

Phys. Rev. Lett. 96, 080403 (2006). 
[18] R. B. Rodrigues, P. A. Maia Neto, A. Lambrecht and S. 

Reynaud, Phys. Rev. Lett. 96, 100402 (2006). 
[19] M. Levin, A. P. McCauley, A. W. Rodriguez, M. T. H. 

Reid and S. G. Johnson, Phys. Rev. Lett. 105, 090403 

(2010). 

[20] K. A. Milton, P. Parashar, N. Pourtolami and I. Brevik, 
Phys. Rev. D 85, 025008 (2012). 

[21] It is worth mentioning however, that even without con- 
sidering bulk objects and new geometries, a repulsive 
force, that scales like 1/r 5 , was also found when the in- 
teraction between an electron in an atomic Rydberg level 
and the ion of the same atom was considered, see e.g., J. 
F. Babb and L. Spruch, Phys. Rev. A 36, 456 (1987). 

[22] J. A. Kong, Electromagnetic Wave Theory, (John Wiley 
and Sons, Inc., 1986). 



